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ABSTRACT 

We examine the effect of magnetic fields on star cluster formation by performing simulations 
following the self-gravitating collapse of a turbulent molecular cloud to form stars in ideal 
MHD. The collapse of the cloud is computed for global mass-to-flux ratios of oo, 20, 10, 5 
and 3, that is using both weak and strong magnetic fields. Whilst even at very low strengths 
the magnetic field is able to significantly influence the star formation process, for magnetic 
fields with plasma (3 < \ the results are substantially different to the hydrodynamic case. In 
these cases we find large-scale magnetically-supported voids imprinted in the cloud structure; 
anisotropic turbulent motions and column density structure aligned with the magnetic field 
lines, both of which have recently been observed in the Taurus molecular cloud. We also find 
strongly suppressed accretion in the magnetised runs, leading to up to a 75% reduction in the 
amount of mass converted into stars over the course of the calculations and a more quiescent 
mode of star formation. There is also some indication that the relative formation efficiency of 
brown dwarfs is lower in the strongly magnetised runs due to the reduction in the importance 
of protostellar ejections. 
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1 INTRODUCTION 



Understanding how stars like the Sun form is one of the key 
questions central to our understanding of the Universe we live 
in. Whilst we have come a long way in this understanding since 
the pioneering work of Jeans, many of the fundamental ques- 
tions such as the rate, distribution and efficiency of star forma- 
tion remain either unknown or are the subject of vigourous debate. 
One of the key areas of uncertainty is whether or not star forma- 
tion is a rapid jlVlac Low & Klessen"2004': 'Hartmann et al.l l200ll : 
lElmegreenl2007h or slow ( Shu et al. 1987; Tan et al. 2006) process, 
central to which is the relative importance of magnetic fields to the 
star formation process. 

Whether or not star formation is rapid or slow, the fact re- 
mains that molecular clouds are observed to contain magnetic fields 
of sufficient strengths that they cannot be ignored in any com- 
plete theory of how stars form from such clouds (e.g. ICrutcheil 
1 19991 : lBoui-keetal.ll200ll : iHeiles & Crut"che3l2005l) . Furthermore 
molec ular clouds ar e observed to contain supersonic turbulent mo- 
tions ( lLarsoiJll98ll) . so the interaction of turbulence and magnetic 
fields is critical to our understanding of the star formation pro- 
cess. This interaction has been the subject of a number of stud- 
ies which have shown that magnetic fields are not effective in 
preventing the rapid dissipation of supersonic turbulence in the 
absen c e of continued driv ing ( Sto ne et alj 1 19981 : iMac Low et al.l 
1 19981 : IVazauez-Semadeni et aL .20o'oir though the presence of 
magnetic fields does change the dynamics of the turbulence 
jPadoan & Nordlund[20o3 : IPadoan et al]|2007h . However, to date. 



there has been only a handful of simulations which have attempted 
to follow the self-gravitating colla pse of a turbu l ent cloud in the 



presence of a magnetic field (e.g. Li et aljr2004l: Nakamural 
l2006l : IVazauez-Semadeni et al j2005l : lTillev & Pudritzll2007h . 



The ability of magnetic fields to provide support against grav- 
itational instability is determined, for an enclosed region of gas 
threaded by a magnetic field, by the ratio of the mass contained 
within the region to the magnetic flux passing through the surface, 
ie. the mass-to-flux ratio, which for a spherical cloud is given by 

^ = ^ (1) 

where M is the mass contained within the cloud volume, "I> is 
the magnetic flux threading the cloud surface at radius R assum- 
ing a uniform magnetic field Bq. The critical value of M/$ be- 
low which a cloud will be supported agains t gra vitational col- 
lapse is given by (e.g.lMouschovias & Spitzerll 19761 : lMestell[l99^ : 



iMac Low & Klesserj|2004l) 



M\ 2ci 



(2) 



where G and fj,o are the gravitational constant and the permeability 
of free space respectively and Ci is a constant determined numeri- 
cally by Mouschovias & Spitzer ( 1976) to be ci ~ 0.53. Star form- 
ing cores with mass-to-flux ratios less than unity are stable against 
collapse ("subcritical") and conversely, cores with mass-to-flux ra- 
tios greater than unity (termed "supercritical") will collapse on the 
free-fall timescale. Throughout this paper we use the mass-to-flux 
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ratio, defined in terms of the critical value, to quantify the mag- 
netic support of a molecular cloud against collapse and we perform 
simulations using initially supercritical clouds (though to varying 
degrees). Giant molecular clouds are generally thought to have 
magnetical ly subcritical envelopes but to b e supercritic a l in th eir 
inner parts fciolek & Mouschoviai h995l) : ICortes et all ^20051) . a 
picture whic h is largely confirmed by observati onal results on both 
large scal es (McKee 1989: 'McKee et alj|l993h and in dense cores 
( ICrutcheAl999 ; Bour ke et al. 2001 ). 

iLi et alj j2004) examined the role of magnetic fields in self- 
gravitating core formation within a turbulent molecular cloud in a 
periodic box. They found that cores formed within a supercritical 
cloud were also locally supercritical by at least an order of magni- 
tude, indicating that a globally supercritical magnetic field does not 
evolve to produce locally magnetically su bcritical cores, contrary 
to some earlier expectations jMestel&Spitzeil 19561) . Furthermore, 
even with supercritical cores they found a central magnetic field 
strength in cores B oc p^^'^ , similar to observations jCrutcheJ 19991) 
which has often been used as an argument for magnetic support in 
molecular cloud cores. They also foun d strong interaction bet ween 
cores and rotationally supported discs. iTillev&PudritibOOTi) also 
performed simulations of self-gravitating collapse in the presence 
of magnetic fields and found that the observed near-critical cores 
could form naturally from a globally highly supercritical cloud, and 
that these cores generally have the same gas-to-magnetic pressure 
ratio, P as the mean /3 in the global cloud. 

Given that magnetic fields always act to oppose gravita- 
tional collapse, it has often been suggested that they may play 
the dominant role in re gulating the star forma tion efficiency 
in molecular clouds (e.g. iKrumholz & Tanl l2007h . though there 
are also several other good can didates for doing so, including 
turbu l ence-inhibited collapse (e.g. Tillev & Pudritzll2004l : iLi et al.l 
l2004lyazauez -Semadeiu et al.ll2005|), feedback from jets and out- 
flows jNakamura & Lill2005l: iLi & Nakamurall2006l) . the dispersal 
of initially unbound clouds " dClark et alj2005l) or radiation feedback 
from the stars themselves. Secondly magnetic fields are often in- 
voked to solve the 'angular momentum problem' in star formation 
vi a the magnetic brakin g of star forming cores. Recent simulations 
bv lPrice & Bate! ( |2007|) have shown that magnetic fields can have 
a dramatic effect on circumstellar disc formation and on fragmen- 
tation to f orm binary systems (these res u lts have since be e n con- 
firmed bv iHennebelle & Tevssieil l2008h . iTillev & Pudrit3 j2007l) 
looked at the effect of the magnetic fields on the star formation 
efficiency in their simulations, though no clear trend was apparent, 
in part due to limitations of the numerical model. 

Despite the apparent importance of magnetic fields to the star 
formation problem, it was therefore somewhat surprisi ng that the 
purely hydrodynamic calculations of iBate et alj ( |2003|) (hereafter 
BBB03) produced largely the 'right answer' in terms of being able 
to reproduce observed the stellar initial mass function (albeit at low 
number statistics) as well as several other observational characteris- 
tics such as the frequency of binary stars and stellar velocity disper- 
sions. In this and subsequent calculations the initial mass function 
is built up due to the competition between dynami cally interacting 
proto stars in order to accrete from the global cloud ( IBonnell & Bate! 
l2006h . Thus, low mass stars are simply those which have been 
quickly ejected f rom multiple systems and thus have on ly a short 
accretion history i lBate et alj20CQl : lBate & BonnelfcOOSi) . whereas 
higher mass stars are those which form and remain at the bottom of 
deep potential wells and b uild up their mass through accretion over 
time teonnell et al.ll 19971) . However subsequent calculations (Bate 



2008, in prep) have established that purely hydrodynamic calcula- 
tions produce an excess of brown dwarfs relative to observations. 

Thus it is crucial to extend these types of calculations, with 
the hope of resolving some of the above issues, by adding the two 
major pieces of missing physics - magnetic fields and the effect of 
radiation transport. This paper presents our first attempt to address 
the former in large-scale simulations. Whilst in a sense the two are 
complementary, since we expect the magnetic field to have an ef- 
fect primarily on larger scales (given the strong physical diffusion 
of the magnetic field on smaller scales), whilst radiation might be 
expected to influence the smaller scale dynamics (ie. fragmenta- 
tion), it is clear that it is imperative to incorporate both pieces of 
physics into these types of calculations. 

The paper presents our first investigation of how magnetic 
fields change the picture of star cluster formation painted by 
BBB03. The numerical method is discussed in ij2]and the initial 
conditions for the simulations are discussed in ^ Results are pre- 
sented in ^and discussed in ^ 



2 NUMERICAL METHOD 
2.1 Hydrodynamics 

We solve the equations of self-gravitating (magneto-) hydrodynam- 
ics using the Smoothed Particle Hydrodynamics (SPH) meth od (for 
recent reviews of SPH see iMonaghanI 12001 iPricd |2004|) . Fluid 
quantities and their derivatives in SPH are evaluated on a set of 
moving particles which follow the fluid motion. The long range 
gravitational force is calculated efficiently using a binary tree algo- 
rithm originally written by Benz et al. ( 1990), although substantial 
modifications have been made to the code since, both in terms of 
efficiency and as improveme nts to the b asic algorithms. Individ- 
ual timesteps were added by iBatd ( |l995l) and sink particles (dis - 
cussed below) were implemented bv Bate. Bonnell & Pried ( Il995h . 
The code at this stage was used for the original BBB03 calcula- 
tions. 

More recently (that is, post BBB03), the hydrodynam- 
ics in the code has been thoroughly updated with state-of-the- 
art SPH algorithms, most notably by adopting the energy and 
entropy-c onservi ng variabl e sm o othing l ength algor it hms d evel- 
oped bv ISpring el & Hern quis^ ^2003) . iMonaghanI ( |2002h and 
IPrice & Monaghan t2004h") and by the i ntroduction of addit ional 
physics in the form of magnetic fields jPrice & Batd |2007|) and 
radiative transfer using the flux-limited diffusion approximation 
dWhitehouse et alj l2005i) (although we do not include radiative 
transfer in this paper). In the variable smoothing length formula- 
tion the density for each particle is calculated according to 



Pi ^^mjWij{hi) 



(3) 



where the smoothing length hi is itself a function of the density in 
the form 



h — ri 



1/3 



(4) 



where 77 is a parameter determining the approximate neighbour 
number (here we choose = 1.2 corresponding to approximately 
60 SPH neighbours). Thus the density summation ^ becomes a 
non-linear e quation for both h and p w hich we solve iteratively as 
described in lPrice & MonaghanI l l2007h . 

Short range gravitational forces (ie. between particles lying 
within each others smoothing spheres) are softened using the SPH 
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kernel with a softening length which is set equal to the SPH 
smoothing length for that particle. We formul ate the force soften- 
ing us ing the formalism presented recently by IPrice & Monaghai] 
( I2OO7I) which ensures that momentum and energy are both con- 
served even though the softening length is a spatially variable quan- 
tity. 



2.2 Magnetohydrodynamics 

The magnetohydrodynamics in the code is based on the rec ent de- 
velopment of MHD in SP H bv IPrice & Monaghad ( l2004jlbh and 
IPrice & Monaghaiil dlOOSt) . For the calculations presented here, as 
in Price & Bate 1 20071) . we use the 'Euler potentials' formulation 



2.3 Equation of state 

The effects of radiative transfer are approximated by adopting an 
equation of state of the form: 



for the magnetic field such that the divergence constraint is satisfied 
by construction. For more details of the Euler potentials formula- 
tion we refer the reader to IPrice & Bate! | |2007|) and to a complete 
descript ion of a code similar (thou gh not identical) to that used here 
given in lRosswog & Pric3 ( l2007h . 

Use of the Euler potentials for the magnetic field evolution is 
a slightly more lim i ted for mulation of MHD than that presented in 
IPrice & Monaghai] ( l2005h in that certain types of initial field ge- 
o metry cannot b e represented in such a formulation (see discussion 
in iRosswog & Price 2007). Whilst this does not present immedi- 
ate difficulties for the simulations presented in this paper (starting 
with a uniform field geometry), we would, for example, not expect 
dynamo processes to be well captured in the Euler potentials for- 
mulation because of the helicity constraint. 

Shocks in both hydrodynamics and MHD are captur ed via dis- 
sipati ve terms corresponding to an artificial viscosi ty (Monaghan 
Il997h and for MHD, resistivity jPrice & Monagharii2004a . .20051) 
with controlling parameters which are individual for each particle 
and evolve with time (thus reducing dissipation away from shocks) 
as described in Price & Monaghan ( 2005) based on the original for- 
mulation of iMorris & Monaghan! ( Il997h . The formulation of arti- 
ficial resistiv i ty in t he E uler potentials' evolution is described in 
IPrice & Bat3 (|2007|) and iRosswog & Pried ( |2007|) . For reference, 
as m IPrice & BatJ ( l2007l) the magneti c force is form ulated using 
the 'Morris formulation' described in IPrice & Mona ghan (2 003) 
which is both stable in the regime where the magnetic pressure ex- 
ceeds the gas pressure whilst conserving momentum sufficiently 
for the accurate simulation of shocks. 

A major limitation to the simulations presented in this prelim- 
inary work is that, at the resolution of the original BBB03 calcula- 
tion (which was determined by the criterion that all of the hydrody- 
namic fragmentation was resolved in the simulation), the artificial 
resistivity plays a dominant role in the evolution of the magnetic 
field on small scales (that is, during the actual collapse to form 
stars). Whilst it may be argued that ideal MHD is also a poor ap- 
proximation for real molecular clouds, it is a wholly undesirable 
situation to have numerical dissipation in place of physical dissipa- 
tion effects. Thus, for example, we are not able at this resolution to 
confidently assert that we have accurately captured the influence of 
th e magnetic field on the fragmentation of individual cores (e.g. as 
m |Price&Batel2007h . Instead we limit ourselves to a discussion of 
the influence of the magnetic field on the large scale structure of the 
cloud and details such as the initial mass function produced in the 
MHD runs should be taken with the appropriate degree of caution. 
Future calculations will be performed at a much higher resolution 
in order to follow the structure of the magnetic field further into the 
collapse. 



P = Kp''. 

where the polytropic exponent 7 is given by 
1, 



7 



P < 10~^^g cm"^ 



7 = 7/5, p > 10""g cm-^ 



(5) 



(6) 



The equation of state is isothermal at low densities (< 10~^^g 
cm~^) where heating and cooling in molecular clouds balance. At 
higher densities the equation of state becomes barytropic with the 
polytropic exponent chosen to match the results of one dimensional 
(spherically symmetric) ca lculations w hich include the full effects 
of radiative transfer ( Masunaga & Inutsuka 2O0Q3 (see BBB03). 

The net effect of the above equation of state is that collapse 
proceeds unhindered until the density reaches the critical density, 
at which point the gas begins to heat as it is compressed, provid- 
ing thermal support which resists collapse. For hydrodynamics this 
critical density sets the mi nimum fra gment mass from which ob- 
jects subsequently accrete ( lBatdl2005l) . 

We caution that the adoption of an equation of state of the form 
l|6j is based on conditions at the centre of a spherically symmetric 
prestellar core of IA/0. Thus this approximation may be expected 
to break down in regions where spherical symmetry is broken - 
most notably this may be true for fragmentation occurring in discs. 
Also, the equation of state in the form ([6f only depends on the lo- 
cal gas density and thus does not account for the propagation of 
radiation which would increase the temperature in the material sur- 
rounding the protostars. Eq uation ^ is also a rather crud e param- 
eterisation even of the ID iMasunaga & Inutsukal ( l200Clh calcula- 
tions, and the effect of changes to the assumed polytropic index has 
not been investigated. The effect of these assumptions on the cloud 
fragmentation is at present uncertain and will ultimately require 
calculations including a self-consistent treatment of radiative trans- 
fer. Fortunately such simulations, whilst not yet including mag- 
netic fields, are starting to be performed (e.g. Iwhitehouse & Batd 
2006) based on the radiative t ransfe r algorithms developed by 
Whitehouse. Bate & Monaghad ^2W^ . 



2.4 Sink Particles 

Sink p articles were introduced into SPH bv lBate. Bonnell & Pried 
l ll995h to enable star formation calculations to be followed beyond 
the point at which stars form in order to study the subsequent global 
cloud dynamics rather than the internal dynamics of the stars them- 
selves. In the calculations described in this paper, sink particles are 
allowed to form (that is, the SPH particle lying closest to the den- 
sity maxima is converted into a sink particle) once the following 
conditions are satisfied: 

(i) the density exceeds 5.5 x 10~^g cm~^ (the exact number 
is somewhat arbitrary - a lower number means that some dynamics 
may be missed whilst a higher number means substantial slowdown 
of the code whilst trying to evolve material at extremely high den- 
sities) 

(ii) a Jeans mass of material is contained within a kernel radius 
(ie. twice the smoothing length) of the particle 

(iii) the material which will form the sink has a ratio of thermal 
to gravitational energy, Ograv < 0.5 and agrav + Pgrav < 1-0 
(where Pgrav is the ratio of rotational energy to the magnitude of 
the gravitational energy). 
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(iv) the divergence of velocity at the particle location is negative 
(ie. material is collapsing) 

(v) all these particles are being moved on the current timestep 

One pitfall in rigourously enforcing the above criteria is that 
it is possible that material within a kernel radius, though part of a 
larger bound object, is not self-bound, unnecessarily blocking sink 
particle formation. A simple method of avoiding this, adopted in 
this paper, is to revert the equation of state to isothermal once a 
density of 10~^^g cm"'' has been reached in order to pass the a 
and P tests and force the formation of a sink particle beyond this 
density. For all practical purposes this is almost identical to simply 
overriding tests ii) and iii) and inserting a sink particle at a density 
of 10~'^'^g cm~^ regardless. It is important to note that this artificial 
change in the equation of state has no effect on the fragmentation 
because the gas never fragments on these scales (ie. in none of the 
calculations do we get sink particles forming close to each other). 

Sink particles, once created, are subsequently allowed accrete 
all material which falls within a fixed accretion radius (set to 5AU 
in the calculations presented here) and is bound to the sink. Short 
range gravitational encounters between sink particles are softened 
using the cubic spline kernel where in these calculations we have 
used a fixed softening length of 4AU for sink-sink interactions 
(note that the softening length for the self-gravitating gas is always 
set equal to the SPH smoothing length and thus varies according to 
the gas density). 

Note that in these calculations , as in previous works 
teate et alj[2003MBate & Bonnellll2005l) , a sink particle is only al- 
lowed to form once the gas has become optically thick to radiation 
and thus no further fragmentation is expected. Thus, provided the 
Jeans mass is sufficiently resolved, the calculations are expected 
to resolve all of the fragmentation present (that is the initial mass 
function is expected to be complete at the low mass end apart from 
extremely hard binaries with orbital separations < 5AU, if indeed 
fragmentation is possible at such scales). 



3 INITIAL CONDITIONS 
3.1 Cloud properties 

The initial molecular cloud is set up identically to that described in 
BBB03: The cloud is initially spherical with a diameter of 0.375 
pc (77,400 AU) and contains a total of 50A/q of material, uni- 
formly distributed, giving an initial density of po = 1.225 x 10~^^g 
cm^^ ("-H2 ~ 3-7 ^ 10^). The cloud free-fall time is t// — 
V37r/(32poG) = 1.90 x 10^ yrs. Thus in observational terms this 
corresponds to a small, relatively dense patch of a molecular cloud 
(although the early evolution whilst the gas remains isothermal is 
scale free). Particles are placed in a uniform random distribution 
cropped to the cloud radius and no particles are placed exterior to 
the cloud, resulting in a significant expansion of the outer layers 
as the collapse proceeds (this would be equivalent to the assump- 
tion of open boundary conditions in a grid-based simulation). The 
cloud is given an initial sound speed of 1.84 x 10"* cm/s, which 
corresponds to a temperature of ~ lOK (the exact temperature de- 
pends on the assumed value for the mean molecular weight). The 
resultant ratio of thermal to gravitational energy is agrav ~ 0.073. 

A total of 3.5 million SPH particles are used in each of the 
calculations (as in BBB03). As previously mentioned, the resolu- 
tion for a cloud of this size was determined in BBB03 by the re- 
quirement that f ragme ntation should be resolved according to the 
lBate&Burker^ ( ll997l) criterion. Given the recent improvements to 



the hydrodynamic algorithm alongside the incorporation of MHD, 
in one sense the purely hydrodynamic calculation presented here 
(see ^ may be viewed as a repeat of the BBB03 calculation with 
a thoroughly updated SPH method, although there are slight differ- 
ences in the initial conditions (discussed below) which mean that 
the hydrodynamic calculation presented here is also not completely 
identical to the original run. 



3.2 Turbulent velocity field 

The cloud is imprinted with a turbulent velocity field as described 
in BBB03 with power spectrum P{k) oc fc^*. Whilst the gener- 
ate d velocity field (p roduced on a grid) is similar to that used in 
the iBate et al] ( l2003h calculations, a slight change in the way in 
which the velocity field was interpolated to the particles has been 
subsequently added to the code and this me ans that, whilst th e hy- 
drodynamic evolution is very similar to the Ib ate et alj ( |2003|) cal- 
culation, it is not identical. The initial velocity field is normalised 
such that the kinetic energy is initially set equal to the gravitational 
potential energy of the cloud, which gives an initial Root Mean 
Square (RMS) Mach number of 6.4 and an initial RMS velocity of 
1.17 X 10^ cm/s. 



3.3 Magnetic fields 

We perform a sequence of calculations with an initially uniform 
magnetic field of progressively increasing strength threading the 
cloud. The strength is parameterised in terms of the mass-to-flux 
ratio expressed in units of the critical value, where we have per- 
formed runs using i\f/"l> — oo (ie. no magnetic field), 20, 10, 5 
and 3. All of the values are supercritical since given the absence 
of ambipolar diffusion in our calculations, subcritical clouds would 
not be expected to collapse. We have verified that this is indeed the 
case by also performing a calculation at a mass-to-flux ratio of unity 
(ie. critical) which, as expected, does not collapse to form stars (the 
cloud flattens along the direction of the magnetic field but under- 
goes a bounce and subsequent expansion rather than collapse). 

Given the cloud dimensions and the mass-to-flux ratio, corre- 
sponding physical magnetic field strengths can be determined for 
each of the runs according to 

where M/<I> is the mass to flux ratio in units of the critical value. 
Thus a run with a critical mass-to-flux ratio would have Bq = 
194/iG and for the runs with mass-to-flux ratios of oo, 20, 10, 5 and 
3 the corresponding field strengths are given by Bo = 0, 9.7, 19, 39 
and 65/10 respectively. 

The magnetic field may also be parametrised in terms of the 
plasma P, the ratio of gas to magnetic pressure, according to 

/3 = 0.0276 \^^) 

(8) 

The five runs presented here thus have initial /3's of oo, 11, 2.8, 0.7 
and 0.25 respectively. Note that the magnetic pressure is dominant 
over gas pressure in the cloud for mass-to-flux ratios < 6 which is 
the case for the two strongest-field runs. Indeed we find that these 
two runs shown far more significant differences compared to the 
weaker field and hydrodynamic runs. 
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Figure 1. Global cloud evolution, shown as column density in the cloud at intervals of 0.2 cloud free-fall times (top to bottom) for the five runs of progressively 
increasing magnetic field strength (left to right), parametrised in terms of the mass-to-flux ratio of the cloud in units of the critical value. Thus Af/$ = oo 
is a hydrodynamic evolution whilst the strongest field run is Af/<I> = 3 (that is, supercritical by a factor of 3). Note the large voids and vertical filamentary 
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The Alfven speed in the initial cloud is given by 

— '»W.(f)-(^)'"(^)""" 

(9) 

giving VA = 0, 7.8 x 10^ 1.6 x 10**, 3.1 x 10* and 5.2 x lO" cm/s 
for the five runs. Thus, the initial turbulent motions in the cloud are 
super-Alfvenic in all cases with Alfvenic Mach numbers of oo, 15, 
7.3, 3.8 and 2.3 respectively. 

The initial magnetic field is defined as a linear gradient in the 
Euler potentials on the particles. Since the g radient of the Euler 
potentials is computed exactly to linear order jPrice & Batell2007l : 
iRosswog & Pricell2007h . the field is thus uniform everywhere (in- 
cluding at the free boundary). As the calculation progresses the 
field is naturally carried by and thus anchored to a surrounding 
medium created by the expansion of the outer layers of the cloud 
(see above). This initial evolution of the field is discussed further in 
i]4.2| and shown for each of the simulations in Figure |4] 

It is worth briefly discussing the validity of starting the cal- 
culation with an initially imposed uniform magnetic field, since 
clearly in reality there will be a mixture of random and or- 
dered components in the field of varying magnitude. However, 
we also start with a uniform density cloud, so density structure 
and non-uniformity in the magnetic field are both generated self- 
consistently by the initially imposed turbulent velocity field (as 
opposed to starting with pre-existing density structure on which 
a magnetic field is imposed). An alternative approach which could 
be explored in future calculations might be to start with a turbulent 
box containing a magnetic field which has been artificially driven 
to a saturated state (in the absence of self-gravity), although even 
in this case it is not clear that this would correspond any better to 
reality, since molecular clouds are clearly not periodic structures 
and the sudden "turning on" of self-gravity is equally question- 
able. Starting with a uniform magnetic field does however provide 
a meaningful upper limit to the effect of the magnetic field on the 
star formation process, since one would expect that any changes in 
the field geometry (for example, using oppositely directed fields in 
different regions or a field with a large random component) would 
tend to decrease the importance of magnetic fields in the star for- 
mation process, as it would be easier for th e fields to reconnect and 
thus dissipate (e.g. lLubow&Pringlelll996l) . 



4 RESULTS 

4.1 Global cloud evolution 

The evolution of the global cloud is presented in Figure [T] which 
shows column density in the cloud at intervals of 0.2 cloud free-fall 
times (top to bottom) for the five different runs in order of increas- 
ing magnetic field strength (left to right). Thus rows correspond to 
snapshots at a fixed time with varying field strength whilst columns 
represent a time sequence at a given field strength. 

The global cloud evolution at early times it/tjf < 0.4) is 
broadly similar in all five cases. Even at t/tff — 1.2 (bottom row 
of Figure [TJ the main distinguishing features of the hydrodynamic 
cloud (overall cloud shape, location of dense regions) remain ap- 
parent down to a mass-to-flux ratio of 10. This is so because for 
mass-to-flux ratios less than 6 (see i|3.3l above), the field does not 
play the dominant role in the gas dynamics of the cloud. However, 
whilst there are striking differences between the hydrodynamic and 
strongly magnetised cases, even in the weaker field runs differences 
due to the magnetic field are apparent. 




Figure 2. Zoomed-in view comparing the outer parts of the cloud in the 
strong field (A//<i> = 3) run (bottom) to the hydrodynamic ran (top) at 0.6 
free-fall times. The strong magnetic field ran shows filamentary structure 
in the column density aligned parallel to the field lines (which are approxi- 
mately vertical - see Figure|4). 



At early times (t/tff < 0.6, top three rows) there are two 
main distinguishing characteristics. The first is that the shock struc- 
ture produced by the initial turbulent velocity field in the dense re- 
gions (which appear yellow in the figure) appears smoother and less 
well-defined than in the hydrodynamic case. We interpret this as be- 
ing due to the additional pressure support given to the cloud by the 
magnetic field. A similar effect is obser ved in hydrodynamic calcu- 
lations when gas pressure is increased ( iBate & Bonn"eill2005 V The 
second notable difference is that the filamentary structure appears 
more filamentary in the less dense regions, particularly evident in 
the strongest field run (Af/cE> = 3) at t/tff = 0.4 - 0.6 (espe- 
cially in the lower parts of the cloud in the figure). This increased 
filamentary structure, or "stripiness" is roughly aligned with the 
large scale magnetic field threading the cloud (see Figure |4j. At 
higher field strengths the field is dominant in these low density re- 
gions of the cloud and thus channels the gas flow along the field 
lines. A close up view of this structure is shown in Figure[2l com- 
paring an enlarged portion of the strongest field run (bottom) at 
t/tff — 0.6 to the same region in the hydrodynamic run (top). A 
similar alignment of filamentary structure with magnetic field di- 
rection has been recently obs erved in maps of the Taurus molecular 
cloud ( [Goldsmith etai]|2005h . 

At later times (t/tff > 0.8) there are further differences in 
the global cloud evolution. The most obvious of these is that in the 
A//$ — 5 and 3 runs large voids are present in the cloud which are 
completely absent from the hydrodynamic calculation (e.g. com- 
paring the rightmost panels of the second last and last rows with 
the hydrodynamic run). These features appear as a result of large 
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scale magnetic flux which remains threaded through the cloud, il- 
lustrated further in Figure |3] which shows a zoomed-in portion of 
the cloud from the A//"I> = 5 run. The plot shows column density 
(top panel) together with a plot of the column-integrated magnetic 
pressure and a map of the integrated magnetic field with strength 
and direction given by the arrows (bottom panel). The single sink 
particle which has formed at this point in this simulation is shown 
in black. Clearly visible is a large void structure to the immediate 
left of the sink, extending to the upper left and diagonally to the 
bottom right in the figure. The lower plot, showing the integrated 
magnetic pressure, appears almost as an inverse of the top panel 
- that is, the column density is low where the magnetic pressure 
is high. Furthermore the magnetic field direction closely traces the 
void structure visible in the column density plot. 

The void is created by material which slides down the mag- 
netic field lines, creating an evacuated region which is unable to be 
refilled by material perpendicular to the field lines. Since the mag- 
netic flux does not change but the gas pressure decreases, the re- 
sult is a region where the magnetic pressure is dominant and which 
prevents material on either side of it from coalescing to form dense 
structures. The magnetic pressure increase in the lower right part 
of Figure [3] is driven further by a "sandwich" compression of two 
dense filaments perpendicular to the field lines (visible above and 
below the void in the top panel), which squeezes the magnetic field 
lines and thus increases the magnetic pressure (essentially until 
the magnetic pressure balances the ram pressure of the filaments), 
whilst increasing the density proportional only to the one dimen- 
sional change in volume. Thus the key to magnetic pressure domi- 
nated void creation is a turbulent velocity field which can produce 
one and two dimensional compressions rather than the isotropic 
compression produced by gravitational forces. 

The ability of the magnetic field to support parts of the cloud 
against collapse has significant implications for the star formation 
rate in the cloud as a whole (and presumably also the overall star 
formation efficiency) and is discussed further in ^4.5\ 

Finally, a delay in the onset and vigour of star formation is 
apparent even in this global view att/tff — 1.2 since stars which 
have been ejected from their parental envelopes are already visible 
in the hydrodynamic cloud at this time whilst none are visible in the 
runs which include a magnetic field. The star formation sequence 
in each case is discussed further in ^43\ below. 

4.2 Magnetic field evolution 

The magnetic field in each of the magnetised runs (A//$ — 
20, 10, 5 and 3) is shown in Figure|4]at intervals of 0.4 cloud free- 
fall times (left to right). In these plots we show streamlines of the 
magnetic field direction (column integrated) in the cloud, overlaid 
on a column-integrated map of the magnetic pressure in the cloud, 
normalised in each case relative to the initial magnetic pressure. 
Thus the colour scale illustrates the relative compression of the field 
in each case. 

In the weaker field runs (A//$ — 20 and 10, top two rows) 
the magnetic field is strongly compressed both by the shocks re- 
sulting from the initial turbulent velocity field (most visible at 
t/tff — 0.4) and subsequently by the gravitational contraction of 
the cloud (t/tff > 0.8, right panels), also resulting in strong dis- 
tortions of the initially straight magnetic field lines. However even 
in the weak field cases, whilst the field is significantly distorted by 
the collapse, the large scale geometry of the field remains imprinted 
into the cloud by the collapse and the net flux threading the cloud 
remains apparent even at late times. In fact the large scale structure 




Figure 3. Close up view of the void structure in the M/<I> = 5 run at 
t/tff = 1.05, showing column density in the cloud (top) and a rendered 
plot of the integrated magnetic pressure with overlaid arrows indicating the 
direction and magnitude of the integrated magnetic field. The lower panel 
is almost an exact inverse of the upper panel, indicating that the magnetic 
field is providing the dominant source of pressure in this region. 



of the field in the outer regions of the cloud is altered very little as 
star formation proceeds in the dense central regions. 

The relative compression of the field decreases as the field 
strength increases (ie. comparing snapshots within the same col- 
umn) and in the stronger field runs (Af/<1> = 5 and 3, bottom two 
rows) the field geometry remains largely uniform as the collapse 
proceeds, with only a relatively small compression of the magnetic 
field. In these cases the magnetic field is able to impart significant 
directionality to the gas motions - particularly in the outer parts of 
the cloud, by channelling material along magnetic field lines (the 
effects of which are clearly visible in the column density plots for 
these runs shown in Figure[T). The anisotropy of turbulent motions 
in the presence of a magnetic field is a clear pre diction of MHD 
turbulence theory (e.g. iGoldreich & Sridhaj|l995l) . This also leads 
to the tantalising possibility that it may be possible to infer and 
thus map both magnetic field strength and direction in molecular 
clouds by measuring anisotropy in interferometric velocity maps 
(Vestuto etal. 2003). 

Finally, Figure |4] also illustrates how the boundary condition 
on the magnetic field is treated in the calculations (discussed above 
in i]3.3t . Initially (left panels) the field is defined only on the par- 
ticles but remains uniform at the boundary because the gradient in 
the Euler potentials is computed exactly to linear order regardless 
of the particle distribution. By t/tff = 0.4 (second panel), the 
outer layers of the cloud have expanded and thus provide an exter- 
nal medium into which the magnetic field remains anchored at later 
times. 
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Figure 4. Evolution of the magnetic field in each of the magnetised runs (M/<^ = 20, 10, 5 and 3, top to bottom) shown at intervals of 0.4 cloud free-fall 
times (left to right). Plots show streamlines of the integrated magnetic field direction overlaid on a colour map of the column-integrated magnetic pressure, 
normalised in each case relative to the initial magnetic pressure (see colour bars). The weaker field runs (top two rows) show strong compression of the 
magnetic field by the gas, whilst in the stronger field cases (bottom two rows) the field is very effective at providing support to the outer regions of the cloud 
where the column density maps show anisotropic structure parallel to the field lines (Figure[T). 



4.3 Star formation sequence 

The star formation sequence in each of the five runs is presented 
in Figures [5] and [6] although it is best appreciated by viewing an- 
imations of each simulatiorQ. The figures show snapshots on a 
zoomed-in portion of the cloud (dimensions 5156 x 5156AU and 
7219 X 7219AU in Figures [5] and [6] respectively) at intervals of 
0.032 free-fall times throughout the evolution. We again caution 
that the results in the MHD cases should be taken as an upper 



limit on the degree of star formation expected in ideal MUD due 
to the high numerical resistivity present on small scales (although 
this may not be completely unrepresentative since we have also ne- 
glected physical diffusion processes such as ambipolar diffusion). 
Nonetheless, the figures serve to starkly illustrate how the effect of 
the magnetic field on large scales can have a significant influence 
on both the degree and manner of star formation which occurs in 
the cloud. 



http://www.astro.ex.ac.uk/people/dprice/research/mcluster/ 
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Figure 5. Close up view of the star formation sequence in eacli of tlie 5 runs, sliown at intervals of 0.032 cloud free-fall times (top to bottom, continued in 
Figure|6) and in order of increasing magnetic field strength (left to right). Whilst the hydrodynamic run collapses in three main regions which subsequently 
merge, the magnetised runs show delayed collapse in some or all of these regions, leading to more quiescent dynamics and an almost complete suppression of 
star formation in the M/ <1? = 3 case. 



4.3.1 Hydrodynamic run 

In the hydrodynamic case (leftmost column), star formation initi- 
ates in three dense cores (two of which are shown in the t/tff = 
1.1 panel, the other collapses to the top left of this figure and is vis- 



ibleatt/t/y = 1.17). The two protostars shown in the f/f// = 1.1 
panel accrete gas rapidly to reach masses of around 0.1 and O.3A/0 . 
The accretion flow forms massive discs around each of these stars, 
both of which subsequently fragment to give a triple and quadruple 
(double binary) system, respectively, from each of which low mass 
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Figure 6. As in Figure|5]but showing the later time evolution (t/t jj > 1.26) in a slightly wider view (dimensions 7219 X 7219AU). The hydrodynamic and 
M/<I> = 20 runs both contain massive disc fragmentation which is delayed in the A//<I> = 10 run and completely absent from the M/<3? = 5 and A//<3? = 3 
runs which show a much more subdued star formation sequence and fewer violent ejections. 



members are ejected (though some "dance" around the combined 
potential of the two systems before being subsequently ejected). 
The two systems (initially separated by ~ 2000 AU) fall towards 
each other on elliptical orbits and eject several protostars at each 
of two periastron passages (seen at closest approach in the panel 
shown at t/tff — 1.13) before the two systems enter a more cir- 



cular orbit around each other (see panel at t/tff = 1.2). The dy- 
namics is then dominated by this "binary" system (that is, with two 
main concentrations of mass orbiting each other) surrounded by a 
circum"binary" accretion flow. One of the stars in one half of the 
binary system briefly forms a circumstellar disc (R ~ 125AU) be- 
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fore it is destroyed by dynamical interactions with low mass stars 
in the process of being ejected from the system. 

The system receives a strong perturbation when the third of the 
first three dense cores (visible in the upper left of the t/tff = 1.2 
panel falling towards the binary system), together with the disc 
which forms around it during the infall, crashes through the main 
system rather like a cannonball and ejects two protostars in a spec- 
tacular "billiard-ball" style encounter. Thus, at t/tff = 1.26 (first 
panel in Figure|6j the dynamics transitions from mostly a two-body 
system (plus perturbations from multiple lower mass members) to 
a "triple" system (ie. three interacting dense cores) surrounded by 
a large and massive circum-triple disc which fragments into mul- 
tiple single and binary systems (subsequent panels in Figure |6] 
t/tff > 1.33). 



4.3.2 M/$ = 20 

At early times (t/tff < 0.6) the highly supercritical run (M/$ — 
20) evolves almost identically to the hydrodynamic case (Fig- 
ure [T] second column), though the filamentary structure appears 
marginally smoother because of the increased pressure provided by 
the (albeit weak) magnetic field. The slight additional pressure pro- 
vided by the magnetic field also changes the initial star formation 
sequence, as two of the three dense cores form slightly later than 
in the hydrodynamic case. More importantly the third core does 
not collapse {t/tff = 1.1 — 1.17 in Figure|5] second column) un- 
til just before it has merged into the main concentration of mass. 
This slight change in the star formation sequence has a dramatic 
effect on the results, since instead of separate multiple systems 
forming (as in the hydrodynamic case), in this case the accretion 
streams coalesce into one very massive disc which subsequently 
fragments (Figure[6l t/tff = 1.26 onwards) and undergoes rapid 
and vigourous star formation (t/tff — 1.26 — 1.39 in Figure[6l(. 
Though we caution that such disc fragmentation may be an artifact 
of the barytropic equation of state employed in the calculations (see 
discussion in i|2.3t , the difference between this run and the hydro- 
dynamic case serves to illustrate the chaotic nature of turbulent star 
formation, in that even the introduction of a weak magnetic field 
can produce a dramatic difference in the results. 



4.3.3 M/$ = 10 

The star formation sequence in the M/<l? — 10 run differs further 
from the hydrodynamic case. At t/tff — 1.17 in this run (cen- 
tral panel of Figure |5} only the densest of the three initial cores in 
the hydrodynamic run has collapsed, though the core visible to the 
upper left of the main star formation at t/tff = 1.17 in the hy- 
drodynamic case can be seen to form, though later (t/tff — 1.23, 
centre panel of last row in Figure|5} and at a much greater distance 
from the densest region. The core which forms to the upper right in 
the hydrodynamic case is completely suppressed in the Af/$ — 10 
run until around t/tff — 1.29 (visible forming in the centre panel 
of the top two rows in Figure [SJ. As a result the subsequent disc 
in the core of the cluster (which fragments spectacularly in the 
M/$ = 20 run) is at earlier times much less massive and forms 
only an accreting binary system (f // — 1.29 in Figure[6}, though 
the mass accretion onto this system at later times (t/tff > 1.33, 
centre panel of last three rows in Figure|6]l causes further fragmen- 
tation. 



4.3.4 M/$ = 5 

The star formation sequence in the A//$ — 5 run (Figure[5l fourth 
column) is almost unrecognisable compared to the hydrodynamic 
case. The first fragmentation occurs in this case in a disc which 
appears edge-on in Figure [5] (t /i// = 1-1 panel) - that is, per- 
pendicular to the global magnetic field direction. Whilst this disc 
fragments to form a multiple system from which a brown dwarf 
is ejected, the subsequent accretion and thus star formation occurs 
at a dramatically reduced level in this run. Thus whilst inside the 
dense core there are still violent interactions between stars, the core 
itself is generally starved of new material by the lower accretion 
rate and is unable to form any more than a handful of stars. Also, 
without the dramatic encounters with protostars formed at greater 
distances, the ejection of low mass protostars is much less efficient. 
An increase in star formation activity occurs ai t/tff > 1.33 as the 
severely delayed collapse of the upper left core occurs (Figure |6l(. 
The role of the magnetic field in suppressing accretion from the 
cloud is quantified and discussed further below. 

4.3.5 M/^ = 3 

In the Af/<1? = 3 run the collapse is strongly channelled along the 
magnetic field lines (Figures [T] and |4]l and there is relatively little 
compression of the global magnetic field by the gas (Figure|4]l. Star 
formation in the cloud (Figure[5] rightmost column) is strongly sup- 
pressed- only a single, isolated core has collapsed by f/i// = 1.23 
(last panel in Figure[5](. It is not until t/tff — 1.29 (Figure[6} that 
a second core collapses and even then the two remain spatially iso- 
lated until around t/tff = 1.39 (last panel of Figure |6ll. The man- 
ner in which star formation proceeds in the cloud is completely 
different from both the hydrodynamic, weak field, and even the 
= 5 case. Features are visible in the column density which 
are clearly not gravitational in origin (for example the "streamer" 
visible at f/f// > 1.29 in Figure [SJ. 

4.4 Initial mass functions 

The initial mass function (IMF) evaluated in each of the runs is 
shown in Figure [T] in order of increasing magnetic field strength 
(left to right, top to bottom), in each case evaluated at 1.5 cloud 
free-fall times. The dark hatched portion of the histogram indicates 
the stars that have ceased accreting from the global cloud whilst the 
light hatched portion indicates those stars tha t are still accret ing at 
t/tff = 1.5. For comparison the slope of the lSalpeteJ jl955 ) IMF 
is plotted (purple line) as well as the IMF s co vering the sub stellar 
population determined byl ChabrieJ ( l2003h and lKroupj JioOIi) . 

Whilst the low number of stars and brown dwarfs formed over- 
all in the simulations precludes a detailed evaluation of the effect of 
magnetic fields on the IMF, there is a hint of some general trends. 
The first is in the overall normalisation - that is, comparing the to- 
tal number of objects formed in each case. For example, comparing 
the strongest field run (M/^ = 3) in Figure |7] (left panel, bottom 
row) to the hydrodynamic and weak field (Af/<I> — 20) cases (top 
row) it is clear that fewer objects are formed overall in the strong 
field case. 

A further question, and one which can only be tentatively 
approached given the limitations of the present simulations, is 
whether or not magnetic fields have an influence on the shape of 
the IMF. The lower right panel in Figure |7] shows the cumulative 
fractional number of stars as a function of mass for each of the 5 
runs of varying field strength. Whilst we caution that the statistics 
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Figure 7. Initial mass functions at 1.5 free fall times for each of the five runs, in order of increasing magnetic field strength (left to right, top to bottom) and 
the cumulative fractional number of stars as a function of mass in all four cases (bottom right panel), with lines corresponding to the hydrodynamic run (black, 
solid), Af/<I> = 20 (red, dotted), A//<I> = 10 (green, dashed), Af/<I> = 5 (blue, long-dashed) and M/^ = 3 (magenta, dot-dashed). The vertical dashed hne 
in each case indicates the characteristic mass in the hydrodynamic run. 



are low, the general direction is a trend towards fewer low mass 
objects with increasing field strength. This is better illustrated by 
looking at cruder statistics such as the ratio of stars to brown dwarfs 
formed in each of the simulations, given in Table [T] From the table 
we see that in the hydrodynamic and weak field runs there tends to 
be an excess in the number of brown dwarfs relative to the number 
of stars, by as much as a factor of 3. In contrast, in the magneti- 
cally dominated runs, we find roughly equal numbers of stars and 
brown dwarfs. We attribute this to the suppression of accretion in 



the stronger field runs, leading to fewer protostars, fewer dynamical 
interactions and thus fewer ejections of low mass objects. 

4.5 Star formation rate 

The effect of the magnetic field in suppressing accretion from the 
global cloud onto the star forming cores is quantified in Figure [8] 
which shows the total mass in stars (that is, the total mass of all gas 
accreted onto sink particles) as a function of time in units of free- 
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Table 1. Ratio of brown dwarfs to stars formed in each of the runs. Whilst 
there are at present only low number statistics for both populations, there is 
some indication of a trend towards fewer brown dwarfs relative to stars in 
the presence of strong magnetic fields. 



fall times for the 5 runs which form stars (ie. up to Af/"I> = 3, with 
runs indicated by the legend). It is clearly apparent from this figure 
that the mass accretion rate strongly anti-correlates with magnetic 
field strength. Even with a very weak magnetic field (A//$ — 20), 
the accretion rate is clearly lower than the hydrodynamic case up 
until the disc fragmentation which occurs at around t/tff = 1.3 
(Figure|5}. 

In the M/$ = 10 case the accretion rate at early times (up to 
^ t/tff = 1.35) is around half of that in the hydrodynamic run - 
M ~ 8.5MQ/tff compared to Mnyd ~ IBMo/f//. The strong- 
field runs (M/"I> = 5 and M/$ — 3) both show very low initial 
accretion rates M ~ 3— 4M0 /t/ / . The difference between the two 
is that the accretion rate in the A//<1? — 5 run increases dramatically 
at around t/tff = 1.25 as two relatively distant regions of the 
cloud undergo gravitational collapse (further out than the regions 
shown in FigureO, whereas this does not occur in the M/$ — 3 
run. In fact the accretion rate between t/tff = 1.25 — 1.46 in 
the M/$ = 5 run, M ~ 14M0/t//, is only slightly lower than 
the average hydrodynamic rate. However the low initial and later 
accretion rates mean that hy t/tff = 1.5 there is around half of the 
mass in stars in the = 5 run compared to the hydrodynamic 

case (~ 4:Mq compared to ~ 8 A/0, or a 50% reduction). In the 
very strong field run the effect is even more dramatic - hy t/tff = 
1.5 there is only around one quarter of the mass in stars compared 
to the hydrodynamic case 2 A/0 compared to ~ 8A^0, or a 75% 
reduction in the mass converted to stars). 

As an illustration of the effect of magnetic fields in prevent- 
ing lower density gas fr om collapsing (and thus t he effect on larger 
scales) as suggested bv lKrumholz & Tarl ilOO'f ) we plot the mass 
above a particular density in the cloud as a function of time in 
Figure [21 The plot shown M{> p) for density values of p = 
10"^^g/cm^ (top panel) and p = 10"^^g/cm^ (bottom panel). A 
similar trend towards lower mass infall rates with increasing mag- 
netic field strength is also visible in these plots. 



5 DISCUSSION 

We have performed a study of how magnetic fields affect the large 
scale collapse of turbulent molecular clouds to form star clusters, 
computing a range of models with mass-to-flux ratios ranging from 
highly to moderately supercritical (with a corresponding range in 
the ratio of gas to magnetic pressure, /3). Whilst even the weakest 
field runs show differences when compared to the hydrodynamic 
case (e.g. lower accretion rates, different star formation sequences), 
strong differences in the gas dynamics were found to be present in 
the runs where the magnetic pressure dominates the gas pressure 
(ie. f3 < 1), in the form of magnetically-supported voids, column 
density striations due to anisotropic turbulent motions and much 








time [free foil times] 

Figure 8. Effect of the magnetic field on the star formation rate in each of 
the 5 runs. The plot shows the total mass accreted onto sink particles as a 
function of time in each of the calculations. A clear trend is visible in which 
the accretion from the cloud is increasingly suppressed as the magnetic field 
strength increases. The strongest field run (A// <& = 3) shows a 75% reduc- 
tion in the total mass accreted onto stars aX t/tff = 1.5 compared to the 
hydrodynamic case. 
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Figure 9. Mass above a particular density in the cloud as a function of free- 
fall time for the 5 runs (lines as indicated in the legend), showing the results 
for density values of 10~^^g/cm^ (top panel) and lO^^^g/cm"^ (bottom 
panel). Again a clear trend is apparent towards lower mass infall rates with 
increasing magnetic field strength. 



lower accretion rates from the global cloud. The implication of 
these results for both observations of star-forming molecular clouds 
and for our theoretical understanding of star cluster formation are 
discussed below. 



5.1 Relevance to observations 

Zeeman measurements of magnetic fields in molecular clouds 
(fCrutcher 1999; Bourke et al. 2001) suggest typical field strengths 
of i? ~ 10/iG (and generally B < 30pG) for regions with 
T ~ lOK, mass-to-flux ratios which are sup ercritical by a factor 
of ~ 2 — 3, marginally (or perhaps firmly, see lPadoan et al.ll2004h 
super-Alfvenic turbulent velocities and P 0.03 — 0.6 (similar 
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values for /3 an d the ratio of turbulent to Alfvenic velocities are 
also inferred bv lHeiles & Trolandliool in the wider Cold Neutral 
Medium). The inference therefore is that the most realistic of our 
calculations are actually the two strongest field runs (A//$ = 3 and 
M/$ = 5). Given that this is the case, we should expect that all 
of the magnetically-driven features observed in these simulations 
to also be present in observed molecular clouds. 

Recent first results from a large-scale mapping project of the 
Taurus molecular cloud complex by Goldsmith et al. (2005) report 
"ring, arc and bubble-like features" and "striated structure... cor- 
related with the magnetic field direction" in ^^CO maps. In fact 
almost the who l e of t he low density material in Taurus in the 
lOoldsmith et alj j2005l) map appears striated parallel to the global 
magnetic field threading the cloud as mapped from polarization 
measurements, as was found to occur in the simulations in the low- 
density outer regions of the clouds. This is suggestive of a low 
plasma j3 in these regions, since we only find striated structure in 
the two strongest field simulations with /3 < 1 (and most promi- 
nently in the Af/$ = 3 case where /3 = 0.3). This is in broad 
agreement with the polarisation measurements of ICrutcheJ fT999l) 
giving lower limits of /3 > 0.06 for regions within Taurus. 

One of the most striking features of the magnetised collapse 
simulations is the appearance of large-scale magnetic-pressure sup- 
ported voids produced by the large scale magnetic flux tubes 
threading the cloud. In these regions the integrated magnetic 
pressure appears to anti-correlate with the column density of 
the cloud (Figure [3t . Turn ing again to the observations of Tau- 
rus, [Goldsmith etaU ( |2005|) report a "very interesting feature" at 
4'' SO"' + 2^^"^ which ' 'appears as a hole", where it appears that 
"some agent has been responsible for dispersing the molecular 
gas". We propose that this may be a magnetic-pressure driven fea- 
ture. The lack of polarization measurements from this region and 
the orientation of the "hole" suggests that the magnetic field should 
be aligned parallel to the line of sight and would therefore be best 
detectable using Faraday rotation measuremen ts. Whilst very few 
such observations exist, the measurements of IWolleben & Reiclj 
( l2004i) give some hint of increased emission in this regions (sug- 
gesting a field which is parallel to the line of sight), though some- 
what ambiguously. Thus further observations are necessary to con- 
firm this picture. 



5.2 Relevance to theory 

One of the primary effects of magnetic fields in the simulations 
is that, even at field strengths which do not prevent collapse, the 
field can have a significant influence on the star formation rate in 
the cloud. For example overall we find that, after 1.5 cloud free- 
fall times, only 4% of the gas has been converted into stars for 
a marginally supercritical collapse (Af/$ = 3) compared with 
16% in the hydrodynamic case. Similar effects of the magnetic 
field on the star formation rate have been found in MHD cal- 
culations of star formation in the presence of driven turbulence 
I Vazquez' s emadeni et al. (2005). In part this can be attributed to 
the simple fact that the magnetic field adds an extra source of pres- 
sure support to the cloud. Thus we would expect that in the present 
context it would be possible to similarly decrease the star formation 
rate by, for example, simply scaling up the turbulent velocity field 
or by increasing the temperature of the cloud. However we would 
also expect the change in the cloud geometry to be rather different 
under these circumstances, since the magnetic field preferentially 
supports material perpendicular to field lines. The degree to which 
this pressure support is anisotropic depends on the ratio of gas-to- 



magnetic pressure, /3 - for example a weak field will be much more 
readily tangled and thus exert a more isotropic pressure, whereas 
a stronger field will exert a pressure that has a much stronger de- 
pendence on the initial field geometry. Furthermore we would not 
expect to find any of the magnetically driven cloud structures, such 
as the filaments in the cloud envelope aligned with the global field 
direction and magnetically supported voids. IPrice & Bate! ( l2007h 
found that the degree to which the magnetic field can be replaced 
by an equivalent increase in thermal pressure is strongly dependent 
on the field geometry. 

The net result of the lower star formation rate in the magne- 
tised runs is also a decrease in the importance of violent interac- 
tions, leading to fewer ejections and therefore also a trend (though 
tentative) towards more massive stars being formed (ie. relatively 
fewer brown dwarfs). Such a trend might have important implica- 
tions for the variation in the IMF in different environments (e.g. 
iKroupallioOll) and as a function of galaxy evolution. In the local 
environment the IMF appears to have a universal shape (at least 
within the observational uncertainties), although this may be be- 
cause there is also remarkable uniformity in the inferred level of 
magnetic field support in observations of local star formation re- 
gions (Crutcher 1999) (namely that, as discussed above, most star 
forming cores appear to be marginally supercritical with super- 
Alfvenic velocity dispersions and (3 ~ 0.1). 



5.3 Limitations and future directions 

Perhaps the most significant limitation of the calculations presented 
in this paper is that, at the resolution employed, we were not able to 
study the effect of magnetic fields on the small scale fragmentation 
of cores, which ma y be important with respect to the formation of 
circumstellar discs dPrice & Batell2007j)^ and in affecting fragmen- 
tation to form binary system s (Mach ida et alj 200^; Price & Batd 
l2007l;lHennebelle & Tevssren200aiMachida et al.l2007l) . A related 
question which would be interesting to examine is how or whether 
the global magnetic field affects the freq uency and/or strength of 
jets and outflows in molecular c louds (e.g. lMatsumoto & Tomisaka 
l2004lBanerjee & Pudritzl2006l) . 

Whilst increasing the numerical resolution will improve the 
resolution of the small scale magnetohydrodynamics, there is not 
strong motivation to do so as the small scale dynamics is also af- 
fected by other physical processes which have not been modelled 
in the present calculations. The most important of these are the ef- 
fects of radiative transfer (ie. replacing the barytropic equation of 
state) and physical diffusion in the magnetic field due to non-ideal 
MHD. 

As mentioned in i]2.3l we expect from preliminary simulations 
which incorporate a self-consistent treatment of radiative transfer 
(in the flux-limited dif fusion approximation) rat her than a bary- 
tropic equation of state dWhitehouse & Batell2006l) that a full treat- 
ment will affect fragmentation where the assumption of spherical 
symmetry is poor and in the region surrounding a collapsing core 
where the temperature can increase because of the propagation of 
radiation from the high temperature central condensation into the 
lower density gas (a process not captured by a barytropic equation 
of state where temperature is proportional to density). These issues 
in particular may be important for fragmentation in discs, which, it 
should be noted is the prime source of much of the star formation 
activity in the hydrodynamic and very weak-field (i\f/$ = 20) 
calculations presented here. It is therefore crucial to quantify the 
degree to which this fragmentation is physical by performing cal- 
culations which explicitly evolve the radiation field. 
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Secondly, whilst flux-freezing is thought to be a good ap- 
proximation for molecular cloud dynamics on large scales, non- 
ideal MHD effects inclu ding ambipolar (ion-neutral) diffusion 
(e.g. iMouschovias & Paleologou 1981 ; Shu et al. 1987 ), the ef- 
fect of finite conductivity and the H all effect are all im portant 
at some level on smaller scales (e.g. IWardle&Nj[T999i) . This 
in particular applies to t he process of core fragmentation (though 
lOishi & Mac Lowll2006l suggest that ambipolar diffusion is unable 
to set a characteristic mass scale in a turbulent flow because of the 
continued propagation of compressive slow MHD waves below the 
ambipolar diffusion scale). It has also been suggested tha t ambipo- 
lar diffusion rates may be enh anced in turbulent flow (Hei tsch et al.l 
|2004 iLi & Nakamurall2004h and therefore that ion-neutral diffu- 
sion may be important also in the earlier stages of collapse. Thus 
it is imperative that the calculations should be extended to in- 
clude non-ideal MHD effects in order to quantify these effects in 
the present calculations (for example Hosking & Whitworth ( 20o3) 
have already implemented a two-fluid scheme for treating ion- 
neutral diffusion in an SPH context which could be used in future 
calculations). Studying physical diffusion processes also requires 
that the numerical diffusion (e.g. due to the artificial resistivity in- 
troduced in order to capture shocks) be reduced to a level below 
the physical diffusion scale, leading to a much more stringent crite- 
rion for fully resolved MHD simulations (of which those presented 
here are not) compared to purely hydrodynamic runs (where it is 
sufficient to resolve the Jean's length in order to capture the frag- 
mentation scale). 
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